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Abstract. The localization properties of electrons moving in a plane perpendicular to a 
spatially-correlated static magnetic field of random amplitude and vanishing mean are inves- 
tigated. We apply the method of level statistics to the eigenvalues and perform a multifractal 
analysis for the eigenstates. From the size and disorder dependence of the variance of the 
nearest neighbor energy spacing distribution, Pw.l{s), a single branch scaling curve is ob- 
tained. Contrary to a recent claim, we find no metal-insulator-transition in the presence of 
diagonal disorder. Instead, as in the uncorrelated random magnetic field case, conventional 
unitary behavior (all states are localized) is observed. 

The eigenstates at the band center, which in the absence of diagonal disorder are believed 
to belong to the chiral unitary symmetry class, are shown to exhibit a /(Q:)-distribution 
for not too weak random fields. The corresponding generalized multifractal dimensions are 
calculated and found to be different from the results known for a QHE-system. 
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1 Introduction 

The localization properties of non-relativistic electrons moving in a two-dimensional 
disordered system in the presence of static random magnetic fields (RMF) with zero 
mean have been discussed controversially for more than a decade. Investigations of An- 
derson localization of 2d electrons with random complex hopping matrix elements date 
back to the early eighties [|| . The actual study of this model was particularly inspired 
by physical situations encountered in high-Tc superconductors and for composite- 
fermions in the quantum Hall effect at half filling factor. 

Now, for lattice models of width M and length L with uncorrelated random mag- 
netic fields the results gradually seem to converge to the notion that all states are 
localized (see, e.g., and papers quoted therein). In the absence of diagonal disor- 
der, the localization length for systems with M odd is believed to diverge at i? = [P, 
where the eigenstates are suggested to become critical |^ . This special state at E = 
is a consequence of the chiral symmetry ^ which can be broken by diagonal disorder 
or, as reported in the case of quasi-ld systems, by application of periodic boundary 
conditions . The random- flux model was shown to exhibit the same symmetry |^ 
as the Gade-Wegner model ||] . The critical behavior is supposed to be similar to that 
of Dirac fermions with random gauge fields 0, ^, |l^ . 
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Recently, a metal-insulator transition has been asserted for the lattice random- 
flux model incorporating spatially- correlated random magnetic fields and additional 
diagonal disorder . This conclusion was based on a finite-size scaling study of the 
longitudinal conductance which was calculated numerically. 

The purpose of the present paper is twofold. First, we check the alleged metal- 
insulator transition by investigating the scaling behavior of the eigenvalue statistics. 
Second, the spatial fluctuations of the eigenstates near energy E — are calculated and 
the scaling of the corresponding moments are analyzed. In the presence of diagonal 
disorder, we flnd only localized states for the spatially-correlated random fiux model. 
Due to the huge localization length the wave functions near the band center show 
multifractal behavior which, however, differs from the quantum Hall case. 

2 Model with correlated random flux 

A model of non-interacting electrons on a square lattice is considered in a one-band 
tight-binding approximation. The spatially-correlated random fluxes are generated as 
in Ref. flltl, 

(j)im)^Aho/llY,RneM-\rm-rn\Vll), (1) 

n 

where ho and Ic are the strength and correlation length of the random magnetic fields. 
The uncorrelated random numbers i?„ are evenly distributed between (—1,-1-1) and 
the Tm, r„ denote the spatial positions of the flux plaquettes considered. The magnetic 
flux per plaquette, 4>{rn), enters the Hamiltonian via the Peierls phase factors in the 
transfer terms from site k to the neighboring sites k' , 

H = Y.^kclc^+ ^exp(i0fefeO4cfe'- (2) 

k (k^k') 

The total fiux through the plaquette labeled by m is given by the sum of the phases 
Okk' associated with its four links. V = 1 \s taken as the unit of energy and the 
lattice constant a = 1 as the unit of length. The diagonal disorder potentials et are 
realized by a set of uncorrelated random numbers taken from the interval (— VT, -l-W^). 
The eigenvalues have been calculated for square lattices of size L = 64, 128, 192. In 
addition, some selected eigenstates were obtained for systems of size up to L = 305. 

3 Results and discussion 

The results of the level statistics are shown in Fig. |^ where the second moment, 
Io{W, L) = 1/2 Pwx{s) ds, of the level spacing distribution, Pw,l{s), is plotted as 
a function of the system size L and disorder strength W. We take the same parameters 
as in [0 so that the amplitude of the correlated random fields ho = 1.0, the correlation 
length Ic = 5.0, and the Fermi energy Ep — —1.0. The eigenvalues were calculated 
within the interval {E-p — 0.25, -Ep -I- 0.25). Many realizations have been considered 
such that the number of accumulated eigenvalues for each pair (L, W) is « 1 • 10^. 
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Fig. 1 Second moment of the level statistics, Io{W, L), versus disorder strength W. The sys- 
tem sizes are L = 64 (O), L — 128 (+), and L = 192 (□). The inset shows the corresponding 
single branch scaling curve. 



We find that for W > 2, Io{W, L) increases with increasing size L. This behavior is 
expected for localized eigenstates which in the limit L — > oo obey the Poisson statistics 
{Ip = 1.0) of uncorrelated energy levels. For < 2 no scale dependence is observed 
for the sizes investigated but Iq{W, L) equals the random matrix result /rmt = 0.59. 
Therefore, the localization length considerably exceeds the system size L. If there 
were a critical point below W — 2 as suggested in [|ll[, Iq{Wc,L) should be scale 
independent at the critical disorder Wc and larger than /rmt- Hence, we have to 
conclude that for the suggested parameters there is no localization-delocalization- 
transition in square systems. The corresponding single branch scaling curve is shown 
in the inset of Fig |l]. This result is similar to the situation found previously for the 



uncorrelated magnetic field case 1 12 



The increase of the conductance for small W in quasi-ld systems reported in pT| ] 
is presumably a finite size effect because for narrow systems the increase in width M 
may raise the conductance due to the additional transport channels. In the limit of 
large M, however, quantum interference will eventually localize the conductance. 

Turning now to the analysis of the eigenstates, we remind that the spatial fluctua- 
tions of critical eigenstates V'£;,l(^) exhibit multifractal properties, i.e., the scaling of 
the appropriate moments, Pq{l,E,L) = X]i(Sreni(0 follow power-laws 

with an infinite set of unrelated exponents T{q) = {q — l)D{q). Alternatively, one 
considers the corresponding /(a)-distribution which for \q\ < 1 can often be approxi- 
mated by /(a) = D{0) — {a — ao)^/[4(ao — D{0))], where in 2d for our lattice system 
D{0) — 2. Hence, ao effectively characterizes the whole distribution. 

In Fig. ||the /(a)-distribution of an eigenstate near £' — of the correlated random 
magnetic fiux model is shown. Similar results are obtained for even L or different Hq 
and Ic parameters as long as the localization length considerably exceeds the system 
size. In most cases, a parabola can be fitted with ao = 2.14 ± 0.03. This number is 
quite distinct from the value ao — 2.29 of a QHE-system [n3|. From the scaling of 
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Fig. 2 The /(Q)-distribution of an electronic eigenstate with energy close to the band 
center for a 2d system with correlated RMF. The parameters are L = 205, W = Q, ho = 1 
and Ic ~ 5. 



the second moment we find a fractal dimension D{2) — 1.87 ± 0.03 for the correlated 
random magnetic field case which also differs from the QHE value 1.6 [ p^ . 

Our values for spatially- correlated RMF are similar to those reported for the uncor- 
related case. The authors of Ref. determined D{2) = 1.8 ±0.06 from the exponent 
5 of the temporal decay of the autocorrelation function using the relation D{2) = 2S 
[ p^ , and D{2) = 1.79 was extrapolated in Ref. from a multifractal analysis similar 
to ours. Multifractal wave functions have also been reported |l^ for massless Dirac 
fermions near E = with a /(Q;)-distribution very close to the one shown here. 

We find two exceptions to the multifractal behavior at E Ki discussed above. 
First, in the presence of diagonal disorder 14^ > 1 a tendency to strong localization 
can be observed which is fully developed for 14^ > 4. Second, if the amplitudes of 
the spatially-correlated {Ic = 5.0) random magnetic fields are weak, hg < 0.2, the 
eigenstates for odd system size L < 305 and W = appear to be extended (ao « 2.04) 
which may be due to the small L or an indication of the chiral unitary symmetry class. 
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